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Abstract 

A new scheme for amplification of coherent gamma rays is proposed. The key elements are 
crystalline undulators — single crystals with periodically bent crystallographic planes exposed to 
a high energy beam of charged particles undergoing channeling inside the crystals. The scheme 
consists of two such crystals separated by a vacuum gap. The beam passes the crystals successively. 
The particles perform undulator motion inside the crystals following the periodic shape of the 
crystallographic planes. Gamma rays passing the crystals parallel to the beam get amplified due 
to interaction with the particles inside the crystals. The term 'gamma klystron' is proposed for the 
scheme because its operational principles are similar to those of the optical klystron. A more simple 
one-crystal scheme is considered as well for the sake of comparison. It is shown that the gamma 
ray amplification in the klystron scheme can be reached at considerably lower particle densities 
than in the one-crystal scheme, provided that the gap between the crystals is sufficiently large. 



PACS numbers: 41.60.-m, 61.85.+P 



I. INTRODUCTION 



Development of coherent radiation sources for a wavelength below 0.1 nm (which cor- 
responds to the photon energy of tens keV or higher), i.e. in the gamma ray range, is a 
challenging goal of modern physics. Such radiation will have many applications in basic 
science, technology and medicine. In particular, they will have a revolutionary impact in 
nuclear and solid state physics as well as in life sciences. 

The present state-of-the-art lasers are capable for emitting electromagnetic radiation from 
the infrared to ultraviolet range of spectrum. X-ray free-electron lasers are currently under 



construction 



|2|]. Moving further, i.e. into gamma ray band, is not possible without new 
approaches and technologies. 

One of the most promising approaches utilizes the spontaneous emission of gamma rays 
by ultrarelativistic charged particles undergoing channeling in periodically bent crystals. 
An undulator-type radiation is emitted in such a system in addition to the well-known 
channeling radiation due to the periodic motion of the particle that follows the shape of 
crystallographic planes or axes. The feasibility of gamma ray generation by ultrarelativistic 
positrons in planar channeling regime in crystals with periodically bent crystallographic 
planes was proven a decade ago {3, 4] (see also the review and references therein). Such 
a device is known as a 'crystalline undulator'. 1 Recently, the feasibility of the crystalline 
undulator utilizing the planar channeling of electrons was demonstrated 

The operation of the crystalline undulator is based on the phenomenon of charged particle 
channeling in crystals 9j (see also the latest review lOj). Channeling takes place if a charged 
particle enters a crystal nearly parallel to a crystallographic axis or plane. It can be confined 
by the electrostatic potential of crystal atoms so that it moves along the axis or plane. It 
is remarkable that the particle follows the shape of the axis or the plane if they are bent. 
Therefore the particle trajectory inside the crystal can be controlled by choosing a suitable 
shape of the axes or planes 

If crystallographic axes and planes are periodically bent, the trajectories of channeling 
particles are approximately sinusoidal. This results into emission of electromagnetic radia- 
tion predominantly in forward direction ^, Q, 0] . The process of emission is very similar to 



1 The term 'crystalline undulator' was introduced in |6| but was not elaborated there. 
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that in the periodic magnetic field of an undulator [l2, Q- However, the electrostatic 
fields inside a crystal are so strong that they are able to steer the particles much more effec- 
tively than even the most advanced superconductive magnets. 2 Due to this fact, the period 
of crystal bending can be made two to four orders of magnitude smaller than the period 
of a conventional undulator. 3 Therefore the particle in a periodically bent crystal will emit 
electromagnetic radiation with much shorter wavelength, i.e. in the gamma ray frequency 
range. 

The emission of the gamma rays will become even more powerful if the positions of the 
charged particles in a crystalline undulator are correlated along the direction of their motion 
in such a way that the waves emitted by all particles have approximately the same phase. 
The resulting radiation is expected to be not only more powerful but also more coherent than 
the spontaneous radiation from uncorrelated particles. The physical picture of the coherent 
emission by correlated beam is similar to what happens in free electron lasers 22], [23, 24. 25]. 
Therefore, we use for this process the term 'lasing effect in a crystalline undulator'. 

First, we consider a one-crystal scheme. The feasibility of obtaining a lasing effect in 

one crystalline undulator has been studied in |4j within the quantum formalism. In the 

current work, this problem is revisited within the classical approach. Indeed, the emission 

of comparatively low energy photons by a high energy particle (i.e., when ftw/e <C 1, with 

huj and e being, respectively, the photon and the particle energy) can be treated classically. 

I [ 

The classical approach [26j has been widely used for the description of conventional free 



electron lasers (see, for instance 



23 



251 ] and references therein). 
Numerical estimation of the parameters of the positron-based crystalline undulator con- 
firm the results of {4], extremely high densities of channeling positrons (~ 10 21 cm -3 
or higher) would be needed to obtain a lasing effect in one periodically bent crystal. This 
suggest the necessity to search for alternative schemes which would be capable of reaching 



2 Indeed, the interplanar electrostatic field in crystals is typically of the order of 10 GeV/cm which is 
equivalent to the magnetic field of approximately 3000 T. The present state-of-the-art superconductive 
magnets produce the magnetic flux density of the order of 10 T with 45 T being currently the highest 
value obtained by combining superconductive and resistive magnets [15j . 

3 Periodical bending can be obtained by the propagation of an acoustic wave along the crystal [3|, ]j| or by 



r-l 



making regularly spaced grooves on the crystal surface either by a diamond blade lg, Il7l ] or by means of 
laser-ablation ]l8(. Another method is a deposition of S13N4 layers onto the surface of a Si crystal 12 1. 
One more possibility is growing of Sii-^Gea; crystals[l9[ with a periodically varying Ge content x 20t 1 2 11] . 
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the lasing effect at considerably lower particle densities. 

Reconsideration of the one-crystal scheme was necessary to develop the classical formalism 
applicable for the description of the lasing effect in crystalline undulator in presence of the 
dechanneling phenomenon (i.e. the process of leaving the channeling regime due to scattering 
by crystal electrons and nuclei). This formalism is more convenient than the quantum one 
for the developing of alternative schemes. Additionally, the one-crystal scheme, being the 
simplest one, provide a good reference point for the comparison with more sophisticated 
set-ups. 

We propose a two-crystal scheme in which the periodically bent crystals separated by a 
vacuum gap. The beam of charged particles passes both crystals successively. In the first 
crystal, a correlation is achieved between the momentum of the particle and its position along 
the beam direction. This correlation is transformed to a density modulation of the beam. 
The gamma radiation is generated mostly in the second crystal by the density-modulated 



beam. This scheme is similar to that of the optical klystron [27|, [28], but is expected to 
operate in the gamma ray range. Therefore we call it a 'gamma klystron'. 

The article is organized as follows. In section [Til a one dimensional model of the one- 
crystal gamma ray amplifier is considered, the formalism is introduced, and numerical es- 
timation are presented. The results of this section are used as a basis for the comparison 
in section [TTTJ where the main subject of our work — the gamma klystron — is considered. 
The results are summarized and discussed in section [V] 

II. ONE-CRYSTAL GAMMA RAY AMPLIFIER 

As it has been already mentioned above, the lasing effect in a crystalline undulator takes 
place if the positions of the channeling particles along the beam direction are correlated in 
such a way that the electromagnetic waves emitted by all particle have approximately the 
same phase. This is accomp 



23 



24j or 'microbunching' 



ished by a spatial modulation (termed usually as 'bunching' 
l|) of the particle density along the direction of the beam 



motion with the period equal to the wavelength of the emitted radiation. To obtain such a 
modulation, initial (seed) radiation from an external source is needed. This radiation may 
be generated by spontaneous emission of charged particles in a crystalline undulator or in 
the field of a infrared laser wave. In both cases the initial radiation has to be well collimated 



4 



I 
I 



L 



i 
i 



FIG. 1: A scheme of the one-crystal amplifier. A charged particle beam (solid lines) and initial 
(seed) radiation (solid wavy lines) enter a crystal with periodically bent crystallographic planes. 
The particle follow the shape of the crystallographic planes and move along nearly sinusoidal 
trajectories (wavy dashed lines). The radiation is amplified due to its interaction with the beam 
in the crystal (see also explanations in the text). 

to ensure sufficient monochromaticity and coherence. 

Under certain conditions which are discussed below, this radiation modulates the density 
of the particles channeling in a periodically bent crystal. Then the bunched beam produces 
additional radiation of the same wavelength. 

In other words, the undulator amplifies the initial radiation. The amplifier is the only 
possible type of lasers in the hard x-ray and gamma ray range. 4 

The scheme of the gamma ray amplifier based on one crystalline undulator (we call it 
'one-crystal gamma ray amplifier') is shown in figure [TJ A charged particles beam and initial 
gamma radiation fall parallel to each other onto a periodically bent crystal of the length L. 
The beam is aligned with the tangent of the bent crystallographic planes at the entrance 
. Therefore, the particles are captured in the channeling mode and move along the planes 
following their shape. In this section we study the conditions under which the interaction 
between the radiation and the channeling particles leads to amplification of the radiation. 

4 In the optical range, the output radiation feedback is commonly used in free electron lasers instead of 
an external source. The feedback is carried out by the mirrors of an optical resonator. Such class of 
devices are called 'oscillators' |25| . This approach cannot be applied in the hard x-ray or gamma-ray 
range because of the absence of the appropriate mirrors. 
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FIG. 2: A positively charged particle trajectory inside a periodically bent crystal. Due to the re- 
pulsive potential of the atomic nuclei (filled circles) the particle gets confined inside the interplanar 
channel of the width d and moves along the bent crystallographic planes following their nearly 
sinusoidal shape with the period A u and the amplitude a and emitting electromagnetic waves. 
Note that the scheme does not reflect the relative scale of d, a and A u adequately. In reality, these 
quantities satisfy the following double inequality d <C a <C A u . 

A. Particle dynamics 

Let us consider a one-dimensional model of the gamma ray amplifier. We assume the 
crystal and the beam to be infinitely wide in the x and y directions, which are perpendicular 
to the direction of the beam propagation (the z axis). 

The crystallographic planes are parallel to the x axes and are bent in ^/-direction, so that 
particles channel in ^-direction and oscillate along the y-axis (see figure [2} ■ 

In what follows the shape y(z) of the periodically bent midplane is assumed to have a 
harmonic form: 

y = a cos(k n z), (1) 
where k u = 2n/\ u with A u being the undulator period. 
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The magnitude of the undulator period exceeds greatly the bending amplitude, i.e. A u ^> 
a (or, equivalently, k u a <C 1) [5|. This inequality, which implies that the crystal deformation 
is an elastic one, follows from the important condition for a stable channeling in a periodically 
bent crystal. Stable channeling occurs if the maximum centrifugal force F c { in the channel 
is less than the maximum interplanar force U / ma:K , i.e. C = F c {/U^ x < 1. For an ultra- 
relativistic particle, the centrifugal force depends on the particle energy as F cf ps e/R m i n , 
where -R m i n = A„/47T 2 a is the minimum curvature radius of the channel with the bending 
profile given by (pQ). Hence, the condition reads [3|, |4| 

C = £^<1. (2) 



A thorough analysis carried out for a number of crystals [5] allowed one to conclude that 
the crystalline undulator is feasible for the values of the bending parameter C lying within 
the interval 0.01 . . . 0.2. For these C values one obtains a ~ 10~ 5 . . . 10~ 3 A u . 

Another important feature of a crystalline undulator concerns the relative magnitude of 
the bending amplitude a and the interplanar spacing d. It was demonstrated (see 0, \i, \^) 
that the operation of a crystalline undulator should be considered in the large-amplitude 
regime, which implies d <C a (typically, a ~ 10 . . . 10 2 <i). This condition, in particular, allows 
one to neglect the oscillations of the particle in the channel. Thus, one can assume that 
under the action of the interplanar field, the particle follows the trajectory ([T]). Formally 
this means that the length s, measured along the trajectory, can be used as a generalized 
coordinate which uniquely characterizes the position of the particle. The corresponding 
conjugate variable is the tangent projection of the particle momentum, 

P> = ~ ^ ( 3 ) 

— s 2 /c 2 



where m is the mass of the projectile and s = ds/dt. 

The evolution of p s is due to the presence of the amplified radiation whose intensity is E. 
The equation of motion reads 

^ = e(E„), (4) 
where e is the charge of the projectile, and n is the unit tangent vector, n = 
(°>2/*> 1 )/\/l + y? (where y' z = dy/dz oc k u a < 1). 

The amplified radiation is sought in the form of a plane wave linearly polarized along the 
y direction and characterized by the wavelength A and the frequency us = ck (here k = 2%/\ 



and c is the speed of light). Neglecting the attenuation 5 of the radiation in the crystal, one 
describes the only non-zero component of the electric field of the amplified wave as follows 

E y = E cos(kz — ujt + 0) . (5) 

Substituting (EJ) into (j4j) and taking into account that ak u <C 1 one derives 
dp s eE Q ak u 



dt 2 



where the following notations are used 



, a 2 k u (k + k u ) . , OI 
sin I ip H sm{2k u xs / 



sin ( ifj + - — — — sin(2fc u xs) — 2k n xs 



(6) 



ip = (k + k u )xs — oot + 0, (7) 

X 



_ (afc u ) 
4 



(8) 



The energy exchange between the particle and the electromagnetic field is most effective 
when the phase tp stays nearly constant, otherwise the first sine term on the right-hand side 
of ([H]) oscillates, thus averaging out the energy exchange 6 . The phase ip is constant provided 
the following resonant condition is fulfilled: 

(k + k u )xs - u = 0. (9) 

Using (jBJ) and expressing the velocity s in terms of the Lorentz factor 7 of the projectile, 
s = c(l — ~ c(l — 7 _2 /2) (the limit 7 3> 1 is assumed), one notices that (jHJ) reduces 

to the following relation between 7 and the wave number k = uj/c: 



k u (4 — a 2 kk u 
Resolving this equation for k one finds 7 



-r=\l , u 2k ^^ ir. (io) 



k ~2 + K^ {U} 



5 It can be shown that the attenuation length L a of gamma rays and hard X rays (see figure 6 in 

is by approximately an order of magnitude larger than the dcchanneling length Ld that determines the 
optimum length of the crystalline undulator. 

6 The phase of the second sine term cannot be made constant. Therefore, this term always oscillates. Hence 
the main contribution to the energy exchange is due to the first sin e te rm. 

7 In accordance with the general theory of free electron lasers (e.g., [23|) the amplification of the electro- 
magnetic wave in an undulator can occur only at the frequencies corresponding to the harmonics of the 
spontaneous undulator radiation, lj^ — kuji (k — 1, 2, . ..). The frequency lo\ of the fundamental har- 
monic coincides with u> = kc which one finds from (jTTJ) . We will consider the emission stimulation in the 
first harmonic only. 



s 



where 



K = 7 r ak u = 27T7 r a/X u (12) 

is the undulator parameter. 

It follows from ffTTj) that k <C k u in the vicinity of the resonance @. Therefore, putting 
k ± & u ~ A; in ([6]), expanding then the right-hand side of the equation in the Fourier series 
and, finally, omitting the oscillating terms one arrives at 

dp s eE ak u . 

— = — J{r])smilj (13) 

where 

r]= Yvh' J(v) = Jo(v)-Ji(ri), (14) 

with Jq,i{v) standing for the Bessel functions. 

Let us derive the equation which describes the evolution of the phase ip (see (JTj)). Dif- 
ferentiating ([3]) one gets dp s /dt = mrfs ~ m7 r 3 s. On the other hand, double differentiation 
of ([7]) yields ip — (k + k u )xs ~ fcs. Combining these equations with (TTH|) and noticing that 
s ~ ct one derives 

^ = -fi 2 sim>), (15) 
which has the form of the pendulum equation. The quantity Q, defined according to 

< 16 » 

has the meaning of the oscillation frequency of the corresponding simple pendulum. 

The derivative dip/ds can be related to the deviation of the particle energy e = , ymc 2 
from its resonance value e T = ^mc 2 : 

# = t = 4fc u 7 - 7, 

In what follows we consider the limit of small gain and small signal. The limit of small 
gain means that the change in the amplitude Eq is much smaller than its initial value 
at the entrance. In other words, the amplitude and, consequently, the frequency Q are 
approximately constant along the undulator. The small signal limit means that E is small 
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enough to ensure the inequality QL <C 1 (L is the length of the undulator). In this case the 
pendulum equation ffToT) can be solved iteratively yielding 

if\„i i /- , fn n2 /sin#o + Cos) sin(^ ) cos(^q)\ 

COO ~ 7- cos(^ + Cos) - cos(^o) + Co (19) 
where ?/> and Co denote the quantities at the undulator entrance. 



B. The gain factor 

The gain factor q(L) characterizes the relative increase in the energy of the amplified 
electromagnetic wave over the undulator length L. This quantity is defined as follows 

«W = ||, (20) 

where £(0) = E 2 /8tt is the radiation energy density at the entrance point, and AS = 
S(L) —£(0) is the increase in the energy density (with £(L) denoting the density at the exit 
from the undulator). 

It follows from the energy conservation law that the radiated energy equals to the decrease 
in the energy of the channeling particles due to the radiative losses. Therefore, to calculate 
the gain factor one can analyze the radiative energy loss of the particles. 

The energy density of the beam particles, £\>{s), at the distance s can be written as 

£ h (s) = (e(s))n(s). (21) 

Here n(s) stands for the volume density of the channeling particles and (e(s)) denotes the 
average energy of a particle at the distance s. The averaging procedure takes into account 
that the initial beam is not spatially modulated, i.e. that for different particles the instants 
of entry into the crystal are not correlated. In other terms, the particles are randomly 
distributed with respect to the phase ip = ip at the entrance point (see (17])). On the 

s=0 

other hand, © suggests that a particle's interaction with the radiation field depends on the 
value of ipo- Therefore, to obtain (e(s)) one averages the energy e(s) = mc 2r y(s) with respect 
to 4>q\ (e(s)) = (27r)" 1 / 27r £:(s)#o- 

In a crystalline undulator the particles move in a medium. Due to collisions with the 
crystal constituents the channeling particle increases its transverse energy e± = c 2 Py/2e. At 
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some point e± exceeds the interplanar potential barrier and the particle leaves the channel. 
The average distance traveled by a particle in a crystal until it dechannels is called the 
dechanneling length In a straight crystal L& depends on the crystal type, the crystal- 
lographic plane, the energy and the type of a projectile. In addition to these, Ld acquires 
the dependence on the parameter C [4| in a bent crystal. The dechanneling effect results 
in a gradual decrease in the volume density of the channeled particles with the penetration 



distance. Roughly, this decrease satisfies the exponential decay law 



29|: 



n(s) = An exp(— s/L d ). (22) 

Here n is the beam density at the entrance of the crystal, and A stands for the channel 
acceptance (i.e., the fraction of the incident particles which is captured into the channeling 
regime) . 

It follows from ( 1211) that the change d£\> of the beam energy density over the interval 
ds contains two terms. One of these, proportional to the derivative of n(s), is due to 
the dechanneling process, whereas another one, proportional to d(e(s)) /ds, describes the 
radiative losses. It is exactly the latter term which is responsible for the change dS of the 
electromagnetic field energy. Therefore, one can write 

d 4- = -mc 2 n(s)^. (23) 
ds ds 

The derivative d^) /ds is calculated using ffT5"]) and f|T7|) . Then, integrating f[2"3"j) . one obtains 

AS = Q 2 [ L n(s) (sin ^(s))ds . (24) 

& Jo 

Using (1T51) and introducing the expansion in powers of fls < QL <C 1, one carries out the 
averaging over if) : 

n 2 ( \ 

(sin(^)) = —2 I sin(Cos) - Co* cos(C s) ) • (25) 
zc, v / 

Substituting (T2"4"l) with fT25l) into (TSUI) one derives the following formula for the gain factor: 
(L,Co) = Stt^tjJ 2 ^) 1 / n(s)(sm(( s) - Co* cos(C s)W (26) 

Here r = e 2 /mc 2 is the classical radius of the particle. The second argument of the gain 
factor in the left-hand side is introduced to indicate its dependence on the quantity £ . 
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In the limit of a short crystal, L C Lj, the dechanneling can be neglected. Then, putting 
n(s) ~ n(0) in (j2Ul) one arrives at 

2 



Co) L« id = ~2nr K^ An V j\ V ) ± (27) 



where u = (qL/2. 



free 



Equation (1271) presents a well-known formula for the gain factor of the conventiona 
electron laser obtained within the small signal and small gain approximation (see e.g. 23|). 

In reality, the dechanneling cannot be neglected. Indeed, the dechanneling length of 
ultra-relativistic positrons (these projectiles are the best candidates for the use in crystalline 
undulators, see discussion in [5|) measured in cm can be roughly estimated as L d ~ O.le 
with the energy in GeV. Hence, the dechanneling length for positrons with e within the GeV 
range does not exceed several millimeters. Therefore, the limit of a long crystal, L ^> is 
of a great interest for an amplifier based on a crystalline undulator. 8 

Assuming in ( 1261) L ^> L d , one extends the upper limit of the integration to infinity and 
carries out the integral. The result reads 

Ld 3 v 
0( L -Co)| L » Ld = lQ7ir k u —An vJ 2 {v) ^ + v2 y > ( 28 ) 

where v = Co-^d- The factor v/(l + v 2 ) 2 attains its maximum of 3y/3/16 at v = l/y/3. 
Therefore, in the limit L ^> L d the maximum gain (with respect to Co) is reached for 
Co = 1/ \/3Ld- Using this value in (I28p one obtains the gain factor g which is optimized with 
respect to Co an d does not depend on the crystal length: 



9 = d(L, Co 



= 3^0^,4^ [ J(v)} 2 ■ (29) 

L»L d 7 r 



Co=l/\/3L d 

This formula coincides up to a numeric factor with the expressions obtained in P, Q within 
the quantum approach. The difference in the numeric factor is due to an additional approx- 
imation made in the cited papers. 



In the consideration we neglect the effect of the radiation attenuation in the crystal. Therefore our results 
are valid if L a L, where L a is the attenuation length. In fact, L a for the photon energies in the range 
10 2 . . . 10 3 keV is on the level of several cm l30j. 
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C. Optimization of the crystal and beam parameters 



The gain factor 0291) depends on a number of parameters. Firstly, it depends on the 
bending period and amplitude. Apart from the obvious dependence of the right-hand side 
on the period via k u = 27rA u , the parameters A u and a enter via rj (see (1121) and fll4H ) . 
Additionally, they enter the acceptance A and dechanneling length L d (as discussed below). 
The latter are also dependent on the type of the crystal and the crystallographic plane as 
well as on the particle energy e r = mc 2r y r . 

In this section we present the scheme which allows one to define the optimum values of 
these parameters, i.e. those which insure the largest possible gain for a given crystal and 
initial beam density n^. Our analysis is restricted to a positron beam. 

To facilitate the analysis it is convenient to express the quantities on the right-hand side 
of (1291) in terms of the parameter C and the undulator parameter K introduced in and 
0T2"]) . respectively. We remind that C stands for the ratio of the maximum centrifugal force 
acting on the particle in the periodically bent channel to the maximum gradient of the 
interplanar potential that keeps the particle in the channel. This definition implies that 
C < 1, otherwise all the particles will be thrown out of the channel by the centrifugal force. 
The lower limit, C = corresponds to a straight channel, a = (or A u = 00). In this limit 
there is no undulator motion and, consequently, there is no radiation amplification. Thus, 
the largest gain corresponds to some optimum value of C lying between and 1. 

The C-dependence of the dechanneling length of an ultra- relativist ic positron can be 



modeled as follows 



L d — (1 — C) 2 L d0 
256 a TF d 

(30) 

9tt z r A 

where L d0 stands for the dechanneling length in the straight channel, otf is the Thomas- 
Fermi radius of the crystal atom, and 

A = ln(^2^mc 2 //) - 23/24 (31) 
rich characterizes the ionization losses of an ultra-relativistic 



is the 'Coulomb logarithm' w' 
particle in amorphous media 



3l| (/ denotes an average ionization potential of the atom). 



Similarly, the acceptance A of a bent channel can be related to the acceptance Aq of the 
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corresponding straight channel 291 ] : 



A — (1 — C)A (32) 



with A 1 — 2a TF /d. 

Expressing the undulator wave number A; u via C and 



r* TP 

k » = l?^f> ( 33 ) 

one re-writes ( 1291) in the form 

g = 3V37rr ^(^) 3 n C(l-Cy\K- 1 V J 2 ( V )] (34) 

mC 2 \ 7 r / L J J? =A' 2 /(4+2if 2 ) 

The formulas flHUj) and (13T1) suggests that the ratio L d0 /j r weakly (logarithmically) de- 
pends on the beam energy: an order of magnitude change in 7 r results in a less than 10% 
change of the ratio. Therefore, we assume this ratio to be constant for a given crystal. Then, 
the gain factor, in addition to the linear proportionality with respect to no, is a function of 
two independent variables, C and K. 

The factor C(l — C) 7 reaches the maximum value of 7 7 /8 8 ~ 0.05 at C = 1/8. The factor 
K^ l r]J'' 2 {ri) 1 as a function of the undulator parameter K, attains its maximum of ~ 0.15 at 
K m 1.2. These are the optimum values of C and K which ensure the maximum gain. The 
latter is given by 

TP T 3 

swx = n) — 5- — r n o- l 35 ) 

mc z 7^ 

The undulator period, A° pt , which corresponds to the optimal values of C and K, one 
derives from ( |33i) : 

>T - 60^. (36) 



max 



The optimal relativistic factor of the beam particles 7 r can be found from ( |T2|) and ( 133!) : 

fT 2 mc 2 mc 2 , , 

j t = = ~ xi. 5 . (37) 

ak u C aU' aU' 



max 



Then the energy Tiw of the emitted photons can be calculated using ( TTTT) : 

2K 3 hmc 3 kmc 3 

fc =™=c(i + jp/2)«»i^ = 18 tf5r- (38) 



max 
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TABLE I: Parameters of gamma ray lasers in the optimum regime for different crystals and planes 
at the temperature T = 4K. The notation a stands for the ratio d/a <C 1. 



Crystal Plane 


d 


U' 

max 


A° pt E 


ftuU 


^0 (flmax 


= 1) 




A 


GeV/cm 


fim GeV 


MeV 


cm" 


3 


C (diamond) (111) 


1.54 


5.16 


59.4 37.7a 132a 2 


1.4 


10 23 


C(graphite) (100) 


3.35 


8.77 


35.0 10.2a 165a 2 


5.3 


10 21 


Si (110) 


1.92 


4.98 


61.6 31.4a 


89a 2 


1.2 


10 23 


Si (111) 


2.35 


6.28 


48.8 20.4a 


47a 2 


4.5 


10 22 


Ge (110) 


2.00 


10.94 


28.0 13.7a 


37a 2 


7.3 


10 22 


Ge (111) 


2.45 


13.55 


22.6 9.1a 


20a 2 


2.9 


10 22 


W (110) 


2.24 


40.74 


7.5 3.3a 


8a 2 


2.0 


10 22 



The wavelength of the produced radiation and the optimal relativistic factor of the beam 
particles 7 r are not fixed by the choice of the optimum values of the parameters K and C. 
They depend on the bending amplitude a. Changing a while keeping the parameters K and 
C constant does not destroy the optimum regime. 

The optimum values of A u and related parameters for different crystals and positron beam 
are shown in Table [B The optimum values of the beam energy e and the photon energy fko 
depend on the ratio a = d/a, which is of the order of 0.1. 

The beam density which is needed to reach g max = 1 has been estimated. As is seen 
from the last column of Table [J extremely high positron densities in the beam are needed 
to obtain a lasing effect in a simple one-crystal amplifier even in optimized regime. This is 
consistent with previous results obtained within the quantum formalism J4]. 

The main reason why an appreciable gain cannot be reached at lower densities is the fact 
that both the beam evolution and the emission of output radiation take place in one crystal 
whose length is limited by a few dechanneling lengths. A deeper density modulation of the 
channeling beam could be obtained in a longer crystalline undulator, but it would not lead 
to the increase of the radiation gain. The reason is the exponential decay of the density of 
the radiating particles because of the dechanneling process. This dilemma is resolved in the 
next section where the scheme of a two-crystal gamma ray amplifier - the gamma klystron 
- is presented. 
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III. TWO-CRYSTAL GAMMA RAY AMPLIFIER THE KLYSTRON 



The scheme of the gamma klystron is shown in figure El Two periodically bent crystals 
of the lengths L\ and L3, respectively, are separated by a vacuum gap of the length L 2 . The 
beam of particles passes the both crystals successively. A correlation between the particle 
momentum and its position along the beam direction is created due to the interaction of the 
channeling particle with the seed radiation in the first crystal. This correlation is further 
transformed into the density modulation of the beam in the vacuum gap. Production of the 
output radiation takes place in the second crystal. 

The idea of the crystal undulator based gamma klystron is very similar to that of optical 
klystron 27|, 



A. Particle dynamics in the Gamma Klystron 

In the first undulator, the phase evolves according to (115p and the solution is given by 
(Fl8|) and (Fl9|) so that the following expressions are valid at exit from the first undulator 



Ci ~ Co + § cosOo + Coil) - t£ cosOq) , 



^1 



Co 



sinQo + CqLi) sinQo) cosQ ) 



(Co^i 



(Coil 



Toil 



(39) 



In vacuum, the particles move along straight trajectories. Therefore, the the equation of 
motion reduces to dPtfj /ds 2 = 0. Hence, 



C2 = Ci , 

at entrance of the second undulator. 

In the second undulator, the phase evolves again according to (fT5l) leading to 

C ~ C2 + 9" COS (^2 + C2S) - ^ C0S(^ 2 
S2 _ t,2 



(40) 



ip ~ 1P2 + C2S + (Sis) 



sin (ifj 2 + (2s) _ sin(ip 2 ) _ cos(-^ 2 ) 
(( 2 s) 2 (( 2 s) 2 C2S 



(41) 



where s is measured from the entrance of the second undulator, < s < L3. 
We assume that 



L 2 3> I/i, I/3. 



(42) 
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FIG. 3: A scheme of the gamma klystron. A beam of charged particles (solid lines) and initial 
(seed) radiation (the solid wavy lines) enter the first crystal with periodically bent crystallographic 
planes of the length L\. The particle follow the shape of the crystallographic planes and move 
along nearly sinusoidal trajectories (the wavy dashed lines) inside the crystal. Interaction between 
the seed radiation and the beam in the first crystal gives rise to a correlation between the particle 
momentum and its position along z axis. This correlation is transformed into a modulation of the 
beam density while the beam travels in the vacuum gap of the length L<i- Then output radiation is 
produced by density-modulated beam in the second crystal of the length L3 (see text for details). 

Therefore, QL2 or Q 2 L\L2 have not to be small even though QLi <C 1 and fls < QL% <C 1. 
Hence, neglecting the terms proportional to QL\ and Qs, but keeping the QL2 terms one 
derives 



sin^ = sin |t/> + Co(^i + L 2 + s) + LiL 2 
Averaging with respect to ipo leads to: 



2 cos(^ + Co^i) - cos(^o) 



(sin ip) 



-Ji 



tt 2 L 2 . {C0L1 
! — - — sm 1 

Co 



sm 



Co^i 



Co ( y + L 2 + s 



(43) 



B. Radiation Gain in the Gamma Klystron 

The increment of the radiation energy density in the second 9 undulator is found as (cf. 

(ED)) 



Ao =!^fjln 2 n{s) (sin ^(s))ds 



k 



(44) 



n 



9 Utilizing of the klystron scheme makes sense only if the radiation gain in second crystal is much larger 
than in the first one. Therefore the gain in the first crystal is neglected. 
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where n(s) stands for the volume density of channeling particles in the second undulator. To 
calculate n(s), one has to take into account the channel acceptance A twice: at the entrances 
of the first and the second undulator 10 . Additionally, the decrease in the number of chan- 
neling particles due to dechanneling in both undulators has to be taken into consideration. 
Therefore, the density of channeling particles in the second undulator reads 

n(s) = A 2 exp ( - — + S ) n . (45) 



La , 

Here no is density of the particles at the entrance of the first undulator. 

We consider the case L3 ^> L& so that the upper limit of the integral in (14411 can be 
replaced by the infinity. One obtains after the integration 

mc 2 7 r 3 „, j9 / L\\ , /^ 9t , sin Mi \ , sin \u\ + 2u 2 + arctan(t>)l . , , 
AS = -^n 2 A 2 n exp -_± U Sl 2 L x L 2 l - L d Li Ui 46 

k \ L d J \ U X J y/l+V 2 

where Ui = ( Li/2, i = 1, 2 and v = (oL d . 

The above expression is derived in the small-gain approximation, but not in the small- 
signal approximation, i.e. it is valid even when the argument of the Bessel function is of 
the order of one. In the following, however, we restrict our consideration to the weak signal 
regime to make a comparison with the one-crystal amplifier. 

In the small-signal case, the radiation gain does not depend on the amplitude of the 
initial wave and has the form 

^£ o_._ LiL 2 L d 2 / Li\ 2 



with 



0(Li,L 2 ,oo) = — = 87rr k u — A n exp -— 77 [J(v)\ h(u 1} u 2 ,v) (47) 

£ 7 r V LdJ 

. sin^! sin \u\ + 2u 2 + arctan(t>)l . . 

h{ui,u 2) v) = — : . (48) 



v 



The optimum value of L\ can be chosen by maximizing the function L\ exp (— Li/L d ). 
The maximum is reached at L\ = L d . 

For fixed values of Li, L 2 and L3, the variables u\, u 2 and v depend only on single 
variable, (q. The optimum value of (0 is found by maximizing the function —h(ui,u 2 ,v). 
Taking into account that u 2 ^> Ui, v due to (|42|) one finds that the maximum is reached at 
u 2 rs — 7r/4 and is approximately equal to 1. 



10 Strictly speaking the acceptances of two undulators are not exactly equal because of different angular 
distributions of particles at the entrances of the undulators. But this difference is small and is neglected 
in our calculations. 
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Therefore the maximum radiation gain which can be reached in the gamma klystron is 
g= (L d ,L 2 ,L 3 )| L 3 «L d =—r k u ^^A 2 n oV [J( V )] 2 . (49) 

( =-n/(2L 2 ) e V 

Here e = 2.7182818 ... is the base of the natural logarithm. 11 

Comparing this formula with ( 1291) one sees that, for the same parameters of the crystals 
and the beam, the radiation gain in the gamma klystron exceeds the gain achievable in the 
one-crystal amplifier by the following factor 

f (klyStr ° n) n = ^U^. (50) 
g (one-crystal) 3V3e L d 

This means that a significant gain can be obtained in the gamma klystron at much lower 
beam densities than in the one-crystal amplifier, provided that L 2 is essentially larger than 



IV. MULTICASCADE AMPLIFIER 



So far we considered a hard X ray or gamma ray amplifier in the small gain regime. 
Estimations in section Hi CI were made for the case g = 1, i.e. when the amplifier doubles 
the intensity of the electromagnetic radiation. Of cause, such a small gain is insufficient. A 
usefull device should be able to amplify the signal by several orders of magnitude. There are 
two methods of increasing the gain. The first one is straightforward. It is merely increasing 
the density of the positron beam. This will switch the amplifier into the large gain regime. 
The potential of this method is, however, limited. For one crystal amplifier, one needs very 
high positron densities, no ~ 10 2 2 cm -3 , to reach even g = 1 . Using the klystron scheme 
will decrease the minimum necessary density by severals orders of magnitude. Nonetheless, 
this quantity will likely remain at the edge of the capabilities of the accelerator technique 
and the sustainability of the crystalline materials. Therefore, increasing the positron density 
much beyond the minimum necessary level will be difficult. 

This difficulty can be overcome by using the second method of increasing the gain: com- 
bining several klystrons into a multicascade amplifier. This scheme is shown in figure HI 
Scattering the beam particles in the crystal decreases the beam intensity and increases its 

11 The constant e in (|49|) and (f50|) should not be confused with the particle charge e. 
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Initial radiation Amplified radiation 




FIG. 4: A scheme of a multicascade hard X ray or gamma ray amplifier. Each cascade consist 
of a klystron fed with a separate positron bunch. The radiation passes the klystron successively. 
If the gain factor of each cascade is g, the radiation intensity will be amplified by the factor of 
G = (1 + g) N , where A" is the number of the cascades. "BD" stands for the positron beam dump. 

emittance. Thereofore, each cascade should be fed by a 'fresh' positron bunch from the ac- 
celerator. The radiation passes each cascade successively and is being gradually amplified. 
The distance between the bunches has to be chosen in such a way that each bunch enters 
the crystal at the same time as the wavepackage of the radiation reaches the corresponding 
cascade. The resulting amplification factor will be G — (1 + g) N , where g is the gain of each 
cascade and N is the number of cascades. 



V. CONCLUSION 



We analyzed two different schemes of the gamma ray amplifier. Both schemes are based 
on crystalline undulators - single crystals with periodically bent crystallographic planes 
exposed to a high energy charged particle beam. Due to the channeling phenomenon, ultra- 
relativistic charged particles move along nearly sinusoidal trajectories inside the crystal. 
Initial gamma radiation traveling through the crystal parallel to the beam is amplified due 
to interaction with the channeling particles. 

We have shown, that the simplest one-crystal scheme would require extremely high (more 
than 10 21 particles per cubic centimeter) positron beam densities to obtain significant am- 
plification. Which makes practical realization of this scheme extremely challenging. 

On the other hand, the two-crystal klystron scheme seems to be more promising. A 
significant gain could be obtained at much lower densities, provided that the distance L 2 
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between two crystals in the klystron exceeds greatly the dechanneling length L d . The ques- 
tion how large the distance L 2 can be in a realistic device cannot be answered within the 
present approach. The simple one-dimensional model considered in the present paper does 
not put any restrictions on the distance L 2 . In reality, the restrictions on L 2 come from 
the energy spread (longitudinal temperature) of the beam. This is, however a technical 
restriction that depends on the quality of the beam determined by the parameters of the 
accelerator. Physical restrictions are more important. They appear, for instance, due to the 
longitudinal velocity spread induced by the channeling oscillations and incoherent scattering 
in the first crystal or the beam divergence in the vacuum gap due to the volume charge. 

Hence further analysis is needed to give the final answer about the feasibility of practical 
realization of the gamma klystron. 

VI. ACKNOWLEDGMENTS 

This work has been supported by the European Commission (the PECU project, Contract 
No. 4916 (NEST)). 



[1] J. Arthur et al. [LCLS Design Study Group Collaboration], SLAC-R-0521 (1998). 

[2] R. Brinkmann, G. Materlik, J. Rossbach, and A. Wagner, DESY-97-048 (1997). 

[3] A.V. Korol, A.V. Solov'yov, W. Greiner, J. Phys. G 24, L45 (1998). 

[4] A.V. Korol, A.V. Solov'yov, W. Greiner, Int. J. Mod. Phys. E 8, 49 (1999). 

[5] A.V. Korol, A.V. Solov'yov, W. Greiner, Int. J. Mod. Phys. E 13, 867 (2004). 

[6] V.V. Kaplin, S.V. Plotnikov and SA. Vorob'ev, Zh. Tekh. Fiz. 50, 1079-1081 (1980), (Sov. 

Phys. - Tech. Phys. 25, 650-651 (1980)). 
[7] M. Tabrizi, A.V. Korol, A.V. Solov'yov, W. Greiner, Phys. Rev. Lett. 98, 164801 (2007). 
[8] M. Tabrizi, A.V. Korol, A.V. Solov'yov, W. Greiner, J. Phys. G 34, 1581 (2007). 
[9] J. Lindhard, Mat.-Fys. Medd. Dan. Vid. Selsk. 34, No. 14 (1965). 
[10] U. I. Uggerh0j, Rev. Mod. Phys. 77, 1131 (2005). 

[11] E. N. Tsyganov, Fermilab preprints TM-682, TM-684 (Batavia, Fermilab 1976) 
[12] V.L. Ginzburg, Izv. Akad. Nauk. SSSR, Ser. Fiz. 11, 165 (1947). 



21 



[13] H. Motz, J. Appl. Phys. 22, 527, ibid. 1217 (1951). 

[14] H. Motz, W. Thon, and R. N. Whitehurst J. Appl. Phys. 24, 826 (1953). 
[15] H.J. Schneider-Muntau, J. Toth, H.W. Weijers, IEEE Trans. Appl. Supercond. 14, 1245 
(2004). 

[16] S. Bellucci et al, Phys. Rev. Lett. 90 (2003) 034801. 

[17] V. Guidi, A. Antonioni, S. Baricordi, F. Logallo, C. Malagu, E. Milan, A. Ronzoni, M. Ste- 
fancich, G. Martinelli, A. Vomiero. Nucl. Inst, and Meth. B 234 (2005) 40. 

[18] P. Balling, K. Kirsebom, D.Q. Le, U.I. Uggerh0j, S. Ballestrero, S. Connell, J. Hartwig, F. 
Masiello. 'Producing and investigating periodic deformations in silicon crystals'. Talk at the 
PECU meeting (Frankfurt am Main, 2007). 

[19] M.B.H. Breese, Nucl. Inst, and Meth. B 132 (1997) 540. 

[20] U. Mikkelsen, E. Uggerhoj, Nucl. Inst, and Meth. B 160 (2000) 435. 

[21] A. V. Korol, W. Krause, A. V. Solov'yov, W. Greiner, Nucl. Inst, and Meth. A 483 (2002) 
455. 

[22] J.M.J. Madey, J. Appl. Phys. 42, 1906 (1971). 

[23] P. Luchini and H. Motz, Undulators and Free-electron Lasers (Oxford University Press, USA 
1990). 

[24] P. Rullhusen, X. Artru and P. Dhez, Novel Radiation Sources Using Relativistic Electrons 

(World Scientific Publishing 1998). 
[25] E.L. Saldin, E.A. Schneidmiller, and M.V. Yurkov, The Physics of Free Electron Lasers 

(Springer, Berlin 2000). 
[26] W.B. Colson, Phys.Lett. A64, 190 (1977). 

[27] N.A. Vinokurov, Proc. 10th Int. Conf. High Energy Charged Particle Accelerators, Serpukhov, 
2 454 (1977). 

[28] N.A. Vinokurov and A.N. Skrinsry, Preprint INP 77-59, Novosibirsk (1977). 

[29] V.M. Biryukov, Yu.A. Chesnokov, V.I. Kotov, Crystal Channeling and its Application at 
High-Energy Accelerators (Springer, Berlin, 1996). 

[30] W.-M. Yao et al, Review of Particle Physics, J. Phys. G 33, 1-1232 (2006). 

[31] V.B. Berestetskii, EM. Lifshitz, L.P. Pitaevski, Quantum Electrodynamics (Pergamon, Ox- 
ford, 1982). 



22 



